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1.  INTRODUCTION 


This  report  documents  a  computational  aerodynamics  and  flight  design  study  of  a  limited 
range  training  round  being  developed  for  U.S.  Army  .50  caliber  ammunition.  The  flight 
requirement  for  the  training  round  is  that  it  must  possess  sufficient  accuracy  for  training 
purposes  to  a  range  of  1000  meters,  with  a  maximum  range  of  2700  meters.  The  requirement 
is  being  satisfied  by  incorporating  a  boattail  design  that  features  sub-caliber  finlets.  The 
finlets  impose  sufficient  roll  damping  on  the  projectile  to  cause  it  to  enter  a  gyroscopically 
unstable,  high-drag  flight  mode  at  some  distance  down  range. 

This  projectile  design  will  satisfy  the  need  for  limited  range  training  ammunition  (LRTA) 
for  the  existing  family  of  .50  caliber  ammunition,  which  includes  (a)  API,  M8  (Armor  Piercing 
Incendiary),  (b)  APIT,  M20  (API  with  Tracer),  (c)  Ball,  M33,  and  (d)  M17  (M33  companion 
tracer  round).  In  addition,  the  LRTA  will  likely  serve  as  a  design  model  for  its  own  companion 
tracer  round. 

The  use  of  a  non-axisymmetric  body  to  impose  spin  control  has  been  reported  (Guidos  & 
Sturek  1987)  for  small  caliber  training  rounds.  That  study  used  unpublished  25-mm  design 
and  free-flight  range  data  provided  by  R.L.  McCoy  of  the  former  U.S.  Army  Ballistic  Research 
Laboratory.  Twisted,  non-axisymmetric  boattail  designs  had  previously  been  examined  and 
later  patented  by  Platou  (1974,  1976)  as  a  means  of  providing  spin  control. 

The  focus  herein  is  mainly  on  the  design  of  the  .50  caliber  LRTA  using  simulation  tech¬ 
niques.  Aerodynamics  coefficients  for  the  LRTA  and  for  existing  .50  caliber  rounds  are 
computed  using  a  three-dimensional  parabolized  Navier-Stokes  (PNS)  viscous  flow  solver. 
The  results  for  the  existing  .50  caliber  rounds  are  validated  against  published  range  data. 
§jx-degree-of-freedom  (6-DOF)  trajectory  simulations  are  made  to  determine  the  yawing 
behavior  and  stable  flight  range  of  the  LRTA,  and  comparisons  are  made  with  Doppler 
radar  test  firing  results.  The  present  LRTA  design  is  evaluated  in  terms  of  the  performance 
requirement,  and  design  modifications  are  recommended. 

2.  CONFIGURATIONS  AND  FLIGHT  CONDITIONS 

Three  existing  U.S.  Army  .50  caliber  projectiles  are  shown  in  Figure  1.  They  are  identified 
as  (a)  API,  M8,  (b)  APIT,  M20,  and  (c)  Ball,  M33.  McCoy  (1990)  has  published  background 
information  about  these  configurations  as  well  as  aerodynamics  results  from  previous  free 
flight  range  tests.  The  range  results  for  the  API,  M8  and  the  Ball,  M33  are  used  in  this  study. 
For  brevity,  the  two  projectiles  are  subsequently  referred  to  as  M8  and  M33.  The  external 
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geometries  of  these  rounds  differ  only  in  minor  surface  details  such  as  rolled  versus  machined 
cannelures.  Figure  2  is  a  schematic  of  the  single  external  geometry  used  to  computationally 
model  both  the  M8  and  M33  projectiles.  The  computational  geometry  does  not  include  any 
of  the  surface  protuberances  found  on  the  actual  rounds. 

The  physical  characteristics  of  the  M33  projectile  are  listed  in  Table  1.  The  M8  and  M33 
projectiles  have  centers  of  gravity  (CG)  locations  that  differ  by  only  about  0.2%  of  their 
total  lengths,  so  differences  in  moment  coefficients  attributable  to  this  CG  difference  are 
negligible. 


Table  1.  Projectile  Physical  Characteristics 


M33 

LRTA 

Reference  Diameter  (mm) 

12.95 

12.95 

Length  (calibers) 

4.46 

5.21 

Mass  (g) 

42.0 

43.7 

CG  (calibers  from  base) 

1.78 

2.70 

/„  (g-cm^) 

7.85 

8.49 

lyy  (g-cm^) 

74.5 

101.8 

The  .50  caliber  LRTA  design  is  shown  schematically  in  Figure  3.  For  simplicity,  the 
nose  section  geometry  is  assumed  to  be  identical  to  the  M33,  although  the  actual  blueprints 
(not  presented  here)  show  subtle  differences  between  the  two  configurations.  The  physical 
characteristics  used  to  model  the  LRTA  design  are  included  in  Table  1.  The  LRTA  boattail 
design  consists  of  a  short  7°  axisymmetric  section  followed  by  a  non-axisymmetric  section  of 
length  slightly  greater  than  one  caliber.  This  non-eixisymmetric  boattail  section  is  formed 
by  engraving  the  cylindrical  section  with  a  mill  cutter  to  form  four  sub-caJiber  finlets. 

Figure  4  shows  the  computational  surface  model  of  the  LRTA  configuration.  For  clar¬ 
ity,  most  of  the  computational  surface  grid  points  eire  omitted  from  the  plot.  The  LRTA 
computational  boattail  model  is  a  simplified  version  of  the  actual  design  in  that  surface 
discontinuities  are  somewhat  smoothed.  This  is  a  characteristic  of  the  fin  geometry  model 
(Rai,  Chaussee,  &  Rizk  1983;  Weinacht,  Guidos,  Sturek  &  Hodes  1986),  used  here  with 
minor  modifications.  The  cross  section  of  the  computational  boattail  geometry  consists  of 
an  axisymmetric  body  which  is  smoothly  coupled  with  spherically  tipped  fins.  In  contrast, 
the  cross  section  of  the  actual  boattcul  section  (refer  to  Figure  3)  is  more  aptly  described  as 
a  square  body  coupled  with  square-tipped  fins.  The  following  approach  was  used  to  provide 
an  appropriate  computational  boattail  geometry  model  for  the  LRTA. 

The  actual  square  boattail  body  is  modeled  using  an  axisymmetric  cross  section  whose 
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(1) 


radius  varies  smoothly  in  accordance  with  the  relation 


nt 


^mtn  +  ^/oc(^re/  ^min) 


in  which  r„/  is  one-half  the  reference  diameter  d,  and  r„..„  is  the  minimum  radius  of  the 
boattail  cross  section.  The  variable  r/„c  is  a  factor  whose  value  decreases  continuously 
and  monotonically  from  one  to  zero.  This  factor  is  defined  using  the  5th  order  polynomial 
expression 


^  fac 


^  10  f  +  6  for  xi^g<x<  x^d  (2) 

V  —  ^htg  )  \  ^end  —  Xbeg  )  \  ®en<f  ^beg  J 


in  which  x  is  the  axial  location  of  the  cross  section  of  interest  (measured  from  the  spherical 
nosetip  vertex),  xy^g  is  the  axial  location  of  the  cylinder-boattail  juncture,  and  Xend  is  the 
axial  location  where  the  boattail  radius  reaches  its  minimum  value.  For  x  >  x^nd,  nt  is  equal 


to  Tmin,  i-e.,  the  cross  section  remains  unchanged. 


Two  different  LRTA  computational  boattail  geometry  configurations  were  generated  by 
assigning  two  sets  of  boattail  geometry  parameters,  and  these  values  are  shown  in  Table 
2.  The  two  configurations  are  intended  to  form  a  geometric  bound  for  the  planform  area 
of  the  finlets  of  the  actual  design.  In  turn,  the  two  configurations  are  expected  to  form  a 
corresponding  bound  on  the  aerodynamic  coefficients  for  the  LRTA.  The  essential  difference 
between  the  two  configurations  is  that  Configuration  #2  has  more  fin  planform  area  than 
Configuration  ^^1. 


Table  2.  LRTA  Boattail  Computational  Geometry  Parameters  (in  Calibers) 


Xbeg 

^end 

^min 

Configuration  #1 

3.75 

4.47 

0.37 

Configuration  #2 

3.75 

4.63 

0.26 

The  service  launch  velocity  for  all  rounds  is  taken  to  be  Mach  2.6.  Atmospheric  sea 
level  free  stream  conditions  are  assumed  for  all  computations.  Computational  results  are 
presented  for  free  stream  Mach  numbers  1.5,  1.96,  and  2.6.  The  Reynolds  numbers  for  these 
velocities  are  taken  as  34.0  million,  44.5  million,  and  59.0  million,  respectively,  based  on  a 
length  of  1  meter.  The  wall  temperature  is  specified  as  294  K,  and  the  flow  is  assumed  to 
be  fully  turbulent. 

Angles  of  attack  are  prescribed  zis  0°  or  2®.  The  roll  orientation  for  the  non-axisymmetric 
configurations  is  such  that  two  finlets  are  aligned  with  the  pitch  plane.  The  magnitudes  of 
the  nondimensional  spinning  and  coning  rates  are  prescribed  in  the  range  from  0.0  to  0.01 
for  various  computational  runs.  The  direction  of  positive  spinning  and  positive  coning  is 
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clockwise  as  viewed  from  the  projectile  base.  The  force,  moment,  and  spin  conventions 
follow  those  put  forth  by  Murphy  (1963). 


3.  PARABOLIZED  NAVIER-STOKES  CFD  TECHNIQUE 


The  PNS  technique  has  been  adapted  and  extensively  used  within  the  Weapons  Technol¬ 
ogy  Directorate  of  the  U.S.  Army  Research  Laboratory  (ARL)  and  is  a  powerful  computa¬ 
tional  fluid  dynamics  (CFD)  research  tool  for  predicting  supersonic  and  hypersonic  projectile 
aerodynamics.  The  PNS  technique  is  a  space-marching  (as  opposed  to  a  time-marching)  tech¬ 
nique;  that  is,  one  numerical  integration  sweep  is  made  from  the  nosetip  of  the  projectile 
to  the  base  to  obtain  a  single  steady  state  solution.  Each  solution  for  the  M33  projectile  at 
angle  of  attack  used  about  5  minutes  of  processing  time  on  a  Cray  X-MP  computer.  Each 
solution  for  the  LRTA  configuration  used  about  15  minutes. 

The  PNS  technology  was  first  applied  to  U.S.  Army  projectiles  to  compute  static  pitch- 
plane  and  Magnus  coefficients  for  spinning  and  non-spinning  shell  and  for  wind  tunnel  models 
at  various  angles  of  attack  (Sturek  &:  Schiff  1981;  Schiff  &  Sturek  1981).  Similar  applications 
were  made  to  shell  at  moderate  angle  of  attack  and  to  finned  KE  projectiles  (Weinacht, 
Guidos,  Kayser,  Sturek  1985;  Weinacht,  Guidos,  Sturek,  Hodes  1986).  The  technique  was 
modified  to  compute  roll  characteristics  of  finned  KE  projectiles  with  exact  fin  geometry 
using  a  rotating  coordinate  frame  (Weinacht  &:  Sturek  1988).  Further  modifications  were 
added  to  compute  the  pitch  damping  of  axisymmetric  (Weinacht,  Sturek,  &  Schiff  1991)  and 
finned  (Weinacht  &  Sturek  1990)  projectiles  using  a  coning  coordinate  frame. 


The  PNS  technique  introduced  by  Schiff  and  Steger  (1979)  is  a  three-dimensional,  finite 
difference,  viscous  flow  solution  procedure  for  attached  supersonic  and  hypersonic  flow  fields. 
The  PNS  technique  spatially  integrates  the  dimensionless,  transformed,  steady,  thin  layer, 
mass-averaged  Navier-Stokes  equations  in  strong  conservation  law  form.  The  governing 
equations  represent  steady  state  conservation  of  mass,  momentum,  and  energy  in  transformed 
coordinates  for  large  Reynolds  number  flows.  The  Cartesian  form  of  the  equations  is 


dE,  OF  dG  ^-idS 


(3) 


These  equations  were  recast  in  cylindrical  coordinates  and  applied  by  Rai  and  Chaussee 
(1983).  The  major  advantage  is  that  the  cylindriceil  coordinate  formulation  requires  only 
three  circumferential  grid  planes  for  axisymmetric  flow  cases  within  the  framework  of  the 
bilateral  symmetry  that  is  imposed.  The  cylindrical  form  of  the  governing  equations  is 
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The  vectors  E,  ,  F,  and  G  contain  the  transformed  inviscid  fluxes.  E,  is  a  modifled  flux 
vector  resulting  from  the  subsonic  sublayer  approximation  (Schiff  ii  Steger  1979).  The  vector 
S  is  the  transformed  vector  of  viscous  terms  that  results  from  the  thin  layer  approximation. 
The  vectors  He  and  Sc  contain  inviscid  and  viscous  source  terms,  respectively,  resulting  from 
the  cylindrical  coordinate  formulation.  The  components  of  the  vectors  for  the  Cartesian 
formulation  are  given  in  many  sources,  including  Schiff  &  Steger.  The  components  of  all  the 
vectors  for  the  cylindrical  formulation  are  given  by  Weinacht  and  Sturek  (1990).  The  three 
transformed  coordinates  are  ^  =  ^(x),  the  axial  (marching);  t]  =  7/(x,  y,  z),  the  circumferential 
coordinate;  and  C  =  ({x,y,z),  the  radial  coordinate.  The  transposed  vector  of  dependent 
variables  is  defined  as 

Q  =  [p,pu,pv,pw,£]  (5) 

in  which  the  density  is  p;  the  axial,  circumferential,  and  radial  velocity  components  are  u, 
V,  and  w,  respectively.  The  total  energy  per  unit  volume  is  £. 

The  solution  is  obtained  at  each  grid  point  using  the  approximately  factored,  implicit, 
delta  form,  flnite  difference  algorithm  of  Beam  and  Warming  (1978).  Second  order  central 
differencing  is  used  in  the  circumferential  and  radial  directions,  and  first  order  one-sided 
differencing  is  used  in  the  marching  direction.  The  solution  is  advanced  downstream  by 
numerically  integrating  in  the  main  flow  direction.  Each  marching  step  requires  a  series 
of  block  tridiagonal  matrix  inversions  (sweeps)  in  the  circumferential  and  radial  directions. 
Fourth  order  explicit  smoothing  terms  are  added  to  suppress  high  frequency  oscillations. 
Second  order  implicit  smoothing  terms  are  added  to  maintain  numerical  stability  in  regions 
of  large  pressure  gradients  (such  as  fin  leading  edges).  Initial  conditions  for  marching  are 
generated  using  the  PNS  method  in  step-back  mode  (Sturek  k  Schiff  1981),  which  assumes 
conical  flow  conditions  near  the  nosetip  and  iteratively  refines  the  solution  to  satisfy  this 
assumption. 

Perfect  gas  behavior  is  assumed.  Turbulence  is  accounted  for  using  the  two-layer,  alge¬ 
braic  eddy  viscosity  model  of  Baldwin  and  Lomax  (1978).  In  that  model,  the  inner  wall 
layer  eddy  viscosity  is  computed  using  a  conventional  Prandtl  mixing  length  with  Van  Dri¬ 
est  damping.  The  outer,  or  wake,  layer  eddy  viscosity  is  based  upon  an  evaluation  of  the 
maximum  moment  of  vorticity  and  its  distance  from  the  wall.  The  calculation  of  the  eddy 
viscosity  is  lagged  by  one  marching  step. 

The  outer  boundary,  which  consists  of  the  bow  shock,  is  shock  fitted  using  the  implicit 


procedure  reported  by  Rai  and  Chaussee  (1983).  The  no-slip  condition  is  enforced  at  the 
body  surface  and  the  pressure  is  held  constant  across  the  subsonic  portion  of  the  boundary 
layer  (i.e.,  the  subsonic  sublayer  approximation).  The  energy  is  defined  from  the  pressure 
using  the  perfect  gas  law.  The  wall  temperature  is  specified,  and  the  density  is  determined 
from  the  temperature  and  pressure. 

4.  RESULTS 

4.1  Aerodynamic  Coefficients  for  M8  and  M33.  Aerodynamics  computations 
were  performed  for  the  M8  and  M33  projectiles  and  the  aerodynamics  coefficients  were 
compared  to  free  flight  range  data  published  by  McCoy  (1990).  The  experimentally  obtained 
aerodynamic  coefficients  were  produced  from  standard  yaw-drag  fits  of  the  equations  of 
free  flight  motion  from  Murphy  (1963).  Since  the  CFD  results  were  generated  at  smedl 
angles  of  attack,  aerodynamic  coefficients  from  range  firings  with  average  angles  of  attack 
greater  than  5°  were  excluded  from  the  comparisons  to  avoid  introducing  large  angle-of- 
attack  nonlinearities  into  the  analysis. 

The  pitch-plane  range  data  are  published  as  the  lift  force  coefficient,  Cl„,  and  the  com¬ 
putational  pitch-plane  results  are  in  the  form  of  normal  force  coeflRcient,  Csa-  The  lift  force 
coefficient  was  converted  here  into  the  normal  force  coefficient  by  employing  the  meaisured 
drag  coefficient,  Cd,  and  the  trigonometric  relationship 

Cn  =  Clcosq  +  Cpsina  (6) 

in  which  Cm  and  Cl  are  normal  and  lift  force  coefficients,  respectively,  at  a  specific  angle  of 
attack.  The  relationship  is  simplified  by  assuming  small  angle  of  attack  into 

Cnc  =  Clc  +  Cd  (7) 

In  Figure  5,  the  comparison  of  pitching  moment  coefficient,  Cm^i  for  the  M8  and  M33 
projectiles  is  shown.  The  computed  value  at  Mach  2.6  is  within  the  scatter  of  the  range  data, 
which  is  about  15%.  The  value  at  Mach  1.5  is  within  about  5%  of  the  data;  the  value  at 
Mach  1.95  is  within  10%  of  the  data.  Figure  6  shows  the  comparison  of  normal  force  center 
of  pressure,  CPm-  The  computed  values  at  Mach  2.6  and  Mach  1.96  are  within  the  scatter 
of  the  data.  The  computed  value  at  Mach  1.5  is  about  0.15  calibers  forward  of  the  measured 
value.  Figure  7  shows  the  comparison  of  normal  force  coefficient,  Csa-  The  computed  values 
are  within  the  scatter  of  the  data  or  within  about  5%  except  at  Mach  1.5,  which  is  off  by 
about  15%. 

Figure  8  shows  the  comparison  of  Magnus  moment  coefficient,  Cmpc-  Most  of  the  range 
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data  indicate  that  the  M8  has  a  positive  Magnus  moment  below  Mach  1.5,  while  the  M33  has 
a  negative  Magnus  moment.  The  computed  values  indicate  a  positive  Magnus  moment  at 
all  three  Mach  numbers.  The  computed  value  at  Mach  2.6  is  within  the  scatter  of  the  data. 
At  Mach  1.96  and  Mach  1.5,  CMpa  is  overpredicted  by  amounts  which  are  about  the  same 
as  the  scatter  of  the  data.  The  range  values  of  Cm^o  below  M=1.5  vary  widely,  indicating 
the  degree  of  difficulty  in  obtaining  repeatable  Magnus  measurements  for  the  M8  and  M33 
projectiles.  Figures  9  and  10  show  the  computed  values  of  Magnus  force  coefficient, 
and  Magnus  force  center  of  pressure,  CPm.  These  coefficients  were  not  obtained  from  the 
range  firings. 

Figure  11  shows  the  computed  and  measured  values  of  pitch  damping  moment  coefficient, 
Cm,  +  Cma-  At  Mach  2.6,  the  computed  value  is  within  the  scatter  of  the  range  data,  while 
at  Mach  1.96  and  Mach  1.5,  the  computed  values  agree  within  about  10%  of  the  data. 

Figure  12  shows  the  computed  values  of  roll  damping  coefficient,  Ci^.  The  roll  damping 
coefficient  was  not  obtained  in  the  range  firings  since  the  configurations  were  not  fitted  with 
roll  pins.  It  is  noted  that  the  roll  damping  is  held  constant  with  respect  to  angle  of  attack 
in  the  6-DOF  simulations  to  be  presented. 

Figure  13  shows  the  computed  and  measured  values  of  zero-yaw  drag  coefficient,  Cdo- 
Three  sets  of  predicted  values  of  Cdq  shown.  The  first  is  the  PNS  result.  The  second 
includes  a  nosetip  contribution  obtained  by  integrating  a  modified  Newtonian  pressure  dis¬ 
tribution  (Anderson  1990).  The  third  includes  a  base  drag  contribution  as  predicted  by  the 
semi-empirical  ‘McDrag’  code  (McCoy  1981).  The  range  results  are  bound  by  the  second 
and  third  sets  of  predicted  drag  values. 

4.2  Aerodynamic  Coefficients  for  LRTA.  Aerodynamics  coefficients  for  the  two 
LRTA  configurations  were  computed  using  the  same  velocities  and  free  stream  conditions  as 
the  M8  and  M33  configurations.  The  LRTA  results  at  Mach  1.5  contain  unwanted  numerical 
oscillations  over  the  final  1/2  caliber  or  so,  and  efforts  to  dampen  the  oscillations  were 
unsuccessful.  These  oscillations  (which  are  shown  in  a  subsequent  figure)  are  not  unexpected 
since  Mach  1.5  is  an  extremely  low  velocity  at  which  to  apply  the  PNS  technique.  The  results 
are  included  here  for  completeness  and  for  future  reference. 

Figures  14,  15,  and  16  show  Cn„  as  a  function  of  x/d  for  Configurations  #1  and  #2 
at  the  three  Mach  numbers  of  interest.  The  computed  results  are  shown  for  the  current 
design  length  of  5.2  calibers  and  for  lengths  extending  approximately  two  calibers  beyond. 
At  the  current  design  length,  Csa  of  Configuration  ^2  is  about  10%  larger  than  that  of 
Configuration  #1.  Figures  17,  18,  and  19  show  CPn  as  a  function  of  xjd  at  the  three  Mach 
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numbers.  At  the  current  design  length,  CPn  of  Configuration  #2  is  about  1/4  caliber  farther 
rearward  than  that  of  Configuration  #1.  The  decrease  of  CyVa  and  CPn  as  boattail  length 
increases  is  a  behavior  that  has  already  been  observed  for  conical  boattail  configurations 
in  supersonic  flow  (Weinacht,  et  al  1985).  CPn,  and  CG  are  input  into  the  6-DOF 

trajectory  simulation  code  and  determine  Cmat  which  influences  the  gyroscopic  stability  of 
the  projectile. 

Figures  20,  21,  and  22  show  Ci^  as  a  function  of  xjd  at  the  three  Mach  numbers.  The 
computational  results  show  the  magnitude  of  C/,,  for  the  LRTA  configurations  to  be  about 
5  times  greater  than  the  M8  and  M33  configurations.  At  the  current  design  length,  Ci^  of 
Configuration  ^2  is  about  30%  greater  than  that  of  Configuration  ^1,  demonstrating  the 
large  sensitivity  of  Ci^  to  fin  planform  area.  The  numerical  oscillations  previously  mentioned 
at  Mach  1.5  are  apparent  in  Figure  22.  In  addition,  the  figures  show  Ci^  to  decrease  with 
respect  to  x/d  over  a  small  portion  of  the  boattail  at  all  three  Mach  numbers  before  increasing 
again.  This  computed  decrease  in  Ci^  with  respect  to  x/d  is  an  unexpected  result.  The 
phenomenon  was  investigated  by  computing  the  flow  field  for  Configuration  ^1  but  with 
two  opposing  fins  removed.  Figure  23  shows  C/p  with  respect  to  x/d  for  the  2-fin  and  4-fin 
configurations  at  Mach  1.96.  The  2-fin  design  does  not  show  a  decrease  in  roll  damping  with 
respect  to  x/d.  It  may  be  concluded  that  computed  decreases  in  C/p  with  respect  to  x/d 
are  attributable  to  the  computed  interaction  of  adjacent  fins  subjected  to  the  expansion  and 
recompression  of  the  flow  over  the  boattail  region. 

Figure  24  shows  the  computed  Cm,  +  Cma  of  the  LRTA  Configurations  #1  and  #2  at 
the  three  Mach  numbers  of  interest.  The  magnitude  of  this  coefficient  is  about  15%  greater 
for  Configuration  ^2  than  for  Configuration  #1  and  is  approximately  twice  that  of  the  M8 
and  M33  projectiles. 

Figures  25  and  26  show  the  computed  Cdq  of  the  LRTA  Configurations  #1  and  ^2, 
respectively,  at  the  three  Mach  numbers  of  interest.  The  predictions  consist  of  the  PNS 
result,  a  modified  Newtonian  nosetip  contribution,  and  a  base  drag  estimate.  The  base  drag 
estimate  was  formed  by  determining  and  integrating  an  average  base  pressure  over  the  base 
of  the  two  configurations.  The  average  base  pressure  of  the  two  LRTA  configurations  weis 
taken  to  be  the  same  as  the  estimated  M33  base  pressure  at  each  Mach  number.  Using  this 
approach,  the  computed  drag  of  both  LRTA  configurations  is  predicted  to  be  less  than  that 
of  the  M8  and  M33  projectiles. 

The  Magnus  force  and  moment  coefficients  are  the  only  aerodynamic  coefficients  not 
presented  for  the  LRTA  but  required  for  the  6-DOF  simulation.  Direct  computation  of  the 
Magnus  coefficients  for  non-axisymmetric  bodies  has  yet  to  be  accomplished  using  CFD 
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technology.  For  typical  non-axisymmetric  projectile  configurations,  such  as  finned  kinetic 
penetrators,  the  Magnus  coefficients  are  not  critical  because  the  spin  rates  are  compara¬ 
tively  low.  The  .50  caliber  LRTA,  however,  uses  spin  rates  large  enough  to  require  consid¬ 
eration  of  Magnus  effects.  A  large  Magnus  moment  could  conceivably  dominate  the  overall 
stability  characteristics  of  the  LRTA,  causing  dynamic  instability  even  though  the  round  is 
gyroscopically  stable. 

The  Magnus  force  and  moment  for  an  axisymmetric  body  is  predominantly  comprised 
of  a  pressure  contribution  that  arises  from  asymmetric  boundary  layer  interaction  with  the 
inviscid  flow  field  (Sturek  k  Schiff  1981).  The  presence  of  fins  gives  rise  to  a  separate 
contribution  in  which  the  force  acts  in  the  opposite  direction  (Platou  1963).  In  the  present 
analysis,  the  computed  Magnus  coefficients  of  the  M33  projectile  are  used  in  the  LRTA  6- 
DOF  trajectory  simulations.  It  is  assumed  that  the  Magnus  contribution  attributable  to  the 
fins  is  zero.  The  Magnus  coefficients  of  the  LRTA  could  probably  be  obtained  from  indoor 
range  firings,  and  a  computational  capability  to  predict  Magnus  moment  coefficient  for  a 
non-axisymmetric  projectile  would  be  valuable  in  a  study  such  as  this. 

4.3  6-DOF  Trajectory  Simulations.  Using  the  predicted  aerodynamic  coefficients, 
6-DOF  simulations  were  performed  using  the  code  whose  development  is  described  by  Fiorellini 
and  Grau  (1992).  The  launch  velocity  is  taken  as  Mach  2.6.  The  quadrant  elevation  is  30°. 
The  initial  roll  rate  is  14627.0  rad/s,  corresponding  to  a  rifled  bore  with  a  twist  of  one 
revolution  per  30  calibers  of  travel.  The  initial  yaw  angle  is  set  to  zero,  and  initial  yawing 
rate  is  set  to  3  rad/s.  Aerodynamic  coefficients  are  supplied  at  the  three  Mach  numbers 
discussed  above,  with  linear  interpolation  between  Mach  numbers  and  linear  extrapolation 
below  Mach  1.5. 

Figure  27  shows  the  computed  Mach  number  versus  range  for  both  LRTA  configurations. 
Both  configurations  exhibit  virtually  the  same  Mach  number  behavior  over  the  first  500 
meters  of  flight.  The  simulation  for  Configuration  #1  reaches  1000  meters  range,  where  the 
Mach  number  is  about  1.5. 

Figure  28  shows  the  computed  gyroscopic  stability  factor,  Sg  (Murphy  1963),  versus  range 
for  both  LRTA  configurations.  The  gyroscopic  stability  factor  is  proportional  to  the  square 
of  the  roll  rate  and  inversely  proportional  to  the  pitching  moment  coefficient.  Following 
Murphy’s  discussion,  a  projectile  is  gyroscopically  stable  if  Sg  is  greater  than  unity.  Con¬ 
figuration  #1  is  launched  with  a  gyroscopic  stability  factor  of  about  1.6  and  decreases  to  a 
value  of  1.0  at  a  range  of  about  750  meters.  Configuration  #2  is  launched  with  a  gyroscopic 
stability  factor  of  about  2.0  and  decreases  to  a  value  of  1.0  at  a  range  of  about  450  meters. 
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The  in-flight  decrecise  of  Sg  for  both  configurations  is  primarily  attributable  to  a  reduction 
in  the  nondimensional  spin  rate,  pd/Voo,  shown  in  Figure  29. 

Figure  30  shows  the  computed  angle  of  attack  versus  range  for  both  LRTA  configurations. 
Both  configurations  axe  computed  to  damp  throughout  the  early  portion  of  flight.  Configu¬ 
ration  ^1  experiences  rapid  yaw  growth  at  a  range  of  about  875  meters,  about  125  meters 
down  range  from  where  it  becomes  gyroscopically  unstable.  Configuration  ^2  experiences 
rapid  yaw  growth  at  a  range  of  about  525  meters,  about  75  meters  downrange  from  where 
it  becomes  gyroscopically  unstable.  Therefore,  the  LRTA  is  predicted  to  damp  throughout 
the  early  portion  of  its  flight  and  experience  rapid  yaw  growth  above  1®  at  a  range  between 
525  and  875  meters.  Results  beyond  this  portion  of  the  flight  were  not  obtainable  since  the 
aerodynamic  coefficients  are  applicable  only  for  small  angles  of  attack. 


4.4  Doppler  Radar  Measurements  &  Discussion.  During  February  1994,  Doppler 
radar  tests  were  conducted  by  the  U.S.  Army  Combat  Systems  Test  Activity  at  Aberdeen 
Proving  Ground,  Maryland  for  the  M33  and  LRTA  configurations.  Doppler  effect  radar 
tracking  system  model  WIOOO  from  Weibel  Inc.  of  Denmark  was  used  to  obtain  the  in-flight 
total  drag  coefficient,  Co-  The  launch  velocity  is  about  Mach  2.6  and  the  quadrant  elevation 
is  25°.  Usable  data  were  obtained  from  approximately  10  rounds  each  for  the  M33  and  LRTA 
configurations. 

Figure  31  shows  a  representative  result  of  measured  in-flight  total  drag  for  the  M33  and 
LRTA.  The  LRTA  experiences  a  rapid  increcise  in  Cp  at  a  range  of  about  1000  meters. 
It  is  assumed  that  this  increase  in  drag,  to  values  an  order  of  magnitude  larger  than  the 
launch  values,  is  caused  by  a  rapid  increase  in  yaw  brought  about  by  the  onset  of  gyroscopic 
instability. 

An  important  observation  of  Figure  31  is  that  the  measured  LRTA  and  M33  drag  coeffi¬ 
cients  are  about  the  same  at  launch.  (These  data  and  the  actual  mean  data  show  the  LRTA 
to  have  slightly  greater  drag  than  the  M33  at  launch.)  The  CFD  predictions  (see  Figures 
13,  23,  and  24)  show  the  opposite  trend  for  zero-yaw  drag,  i.e.,  that  the  M33  has  a  higher 
zero-yaw  drag  than  the  LRTA.  Since  the  yaw  drag  is  included  in  the  Doppler  measurement,  a 
direct  comparison  of  drag  cannot  be  made  between  the  measured  and  predicted  drag  values. 
However,  if  the  measured  and  predicted  drag  trends  are  both  correct,  then  it  may  be  con¬ 
cluded  that  the  LRTA  experiences  larger  initial  launch  disturbances  than  the  M33.  These 
disturbances  would  result  in  larger  initial  yaw  levels  and  increased  total  drag  at  launch. 

Figure  31  also  shows  that  the  LRTA  drag  coefficient  grows  significantly  even  before  it 
reaches  1000  meters  range.  The  growth  is  noticeably  larger  than  that  of  the  M33,  suggesting 
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that  the  LRTA  is  experiencing  yaw  drag  effects  not  predicted  in  the  6-DOF  simulation.  Two 
possible  explanations  are  as  follows:  (1)  the  yaw  levels  in  the  actual  round  are  not  damping 
as  fast  as  those  in  the  6-DOF  simulation  and  (2)  the  yaw  levels  are  actually  growing  (i.e., 
the  round  is  dynamically  unstable). 

A  possible  source  of  dynamic  instability  in  the  LRTA  is  the  Magnus  moment.  As  previ¬ 
ously  mentioned,  direct  prediction  of  the  Magnus  moment  for  the  LRTA  was  not  performed. 
Instead,  the  predicted  Magnus  moment  for  the  M33  was  used  in  the  6-DOF  simulation.  The 
sensitivity  of  the  dynamic  stability  to  Magnus  effects  was  briefly  studied  here.  The  projectile 
is  dynamically  stable  if  the  following  relationship  between  the  dynamic  stability  factor,  Sd, 
and  the  gyroscopic  stability  factor  is  satisfied  (Murphy  1963): 

S42  -  Si)  >  i  (8) 

The  dynamic  stability  factor  depends  primarily  on  the  Magnus  moment  and  pitch  damp¬ 
ing  moment  coefficients.  The  sensitivity  of  the  dynamic  stability  to  the  Magnus  moment 
coefficient  was  studied  by  examining  a  stability  coefficient,  Sd”,  which  reflects  the  relation¬ 
ship  between  the  two  stability  factors.  The  stability  coefficient  is  defined  here  as 

Sj-  =  Sj(2-Sj)-i  (9) 

where  the  projectile  is  dynamically  stable  if  Sd”  is  positive,  in  which  case,  the  maximum  yaw 
is  expected  to  decrease  with  respect  to  range. 

The  stability  coefficient  was  calculated  for  a  parametric  variation  of  Magnus  moment 
coefficient  at  ranges  of  0  and  560  meters  (corresponding  to  in-flight  Mach  numbers  2.6  and 
1.96,  respectively).  The  in-flight  velocities,  spin  rates,  and  aerodynamic  coefficients  from 
the  6-DOF  simulation  (Configuration  #1)  were  used.  Figure  32  shows  that  the  Magnus 
values  used  in  the  6-DOF  simulation  provide  dynamic  stability.  However,  the  figure  also 
shows  that  possible  fin  Magnus  moment  contributions  (which  will  be  negative  in  sign)  of 
equal  magnitude  to  the  Magnus  moment  values  used  in  the  6-DOF  simulation  could  cause 
dynamic  instability. 

It  cannot  be  conclusively  determined  from  the  radar  tests  alone  whether  the  current 
LRTA  design  possesses  sufficient  accuracy  for  training  purposes  at  a  range  of  1000  meters. 
However,  it  is  noted  that  recent  dispersion  test  results  (not  included  here)  show  that  (1)  the 
M33  forms  a  tighter  impact  group  than  the  LRTA  (2)  some  of  the  LRTA  projectiles  show 
considerable  yaw  at  1000  meters,  and  (3)  the  mean  point  of  impact  is  within  1  mil. 

The  aerodynamic  coefficients  presented  here  provide  sufficient  guidance  for  possible  de¬ 
sign  modifications.  One  recommended  modification  is  shortening  the  boattail,  which  would 
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reduce  the  roll  damping  and  increase  the  stable  flight  range  of  the  projectile.  A  shorter 
boattail  would  also  reduce  the  fin  Magnus  effect  and  would  likely  reduce  initial  launch  dis¬ 
turbances,  resulting  in  improved  dispersion  characteristics  over  the  present  design. 

Finally,  it  is  noted  that  the  Doppler  radar  results  showed  the  maximum  range  of  the 
LRTA  to  be  about  1800  meters,  well  within  the  2700  meter  maximum  range  requirement. 
For  comparison,  it  is  noted  that  the  radar  showed  the  maximum  range  of  the  M33  to  be 
about  6000  meters. 


5.  CONCLUSION 

A  computational  aerodynamics  and  flight  design  study  for  .50  Ccdiber  limited  range  train¬ 
ing  ammunition  (LRTA)  has  been  presented.  The  LRTA  possesses  sub-caliber  finlets  that 
impose  sufficient  roll  damping  on  the  projectile  to  cause  it  to  enter  a  gyroscopically  unsta¬ 
ble,  high  drag  flight  mode  at  some  distance  down  range.  A  three-dimensional  parabolized 
Navier-Stokes  (PNS)  viscous  flow  solver  was  applied  to  compute  the  aerodynamic  coefficients 
for  the  LRTA  and  for  existing  .50  caliber  M8  and  M33  projectiles.  The  PNS  results  were 
validated  against  published  range  data.  The  computed  aerodynamic  coefficients  were  used 
in  six-degree-of-freedom  (6-DOF)  trajectory  simulations  to  determine  the  yawing  behavior 
and  stable  flight  range  of  the  LRTA.  Doppler  radar  test  firing  measurements  for  the  LRTA 
were  presented. 

The  PNS  predictions  for  the  M8  and  M33  configurations  agree  with  the  range  data 
mostly  within  the  experimental  scatter.  PNS  predictions  for  the  LRTA  were  made  for  two 
computational  geometries  expected  to  bound  the  aerodynamic  behavior  of  the  actual  LRTA, 
which  consists  of  a  relatively  complex  and  discontinuous  surface  geometry.  The  roll  damping 
coefficient  of  the  two  configurations  is  predicted  to  differ  by  about  30%,  demonstrating  the 
large  sensitivity  of  this  coefficient  to  small  changes  in  fin  planform  area  for  the  LRTA. 

The  6-DOF  trajectory  simulations  show  significant  differences  in  the  stable  flight  range 
for  the  two  computational  configurations.  Configuration  #1  is  predicted  to  maintain  low 
angle-of-attack  flight  to  a  range  of  about  875  meters,  while  Configuration  ^2  is  predicted 
to  maintain  low  angle-of-attack  flight  to  a  range  of  about  525  meters.  The  Doppler  radar 
measurements  show  the  actual  LRTA  to  maintain  low  angle-of-attack  flight  to  a  range  of 
about  1000  meters.  The  Doppler  radar  measurements  also  show  evidence  that  the  LRTA 
experiences  larger  yaw  levels  than  the  M33,  as  well  as  possible  dynamic  instability  because 
of  Magnus  effects  in  the  first  1000  meters  of  flight.  The  measurements  show  the  maximum 
range  of  the  current  LRTA  design  to  be  about  1800  meters,  well  within  the  requirement  of 
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2700  meters. 

The  aerodynamic  coefficients  presented  here  provide  sufficient  guidance  for  possible  de¬ 
sign  modifications.  One  recommended  modification  is  a  shortening  of  the  boattail,  which 
would  reduce  the  roll  damping  and  extend  the  stable  fiight  range  of  the  projectile.  A  shorter 
boattail  would  also  reduce  the  fin  Magnus  effect  and  would  likely  reduce  initial  launch  dis¬ 
turbances,  resulting  in  improved  dispersion  characteristics  over  the  present  design. 

The  use  of  sub-caliber  finlets  for  .50  caliber  limited  range  training  ammunition  appears  to 
be  a  viable  concept.  The  present  .50  caliber  LRTA  design  is  close  to  satisfying  the  flight  re¬ 
quirements  for  this  mission.  Such  a  non-axisymmetric  boattail  concept  is  also  being  pursued 
for  application  to  5.56-mm  limited  range  training  ammunition.  Overall,  the  flight  design 
analysis  presented  in  this  study  demonstrates  how  state-of-the-art  CFD  technology  and  ex¬ 
perimental  development  programs  together  form  a  powerful  projectile  design  capability. 
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Figure  1.  Photograph  of  (a)  API,  M8,  (b)  APIT,  M20,  and  (c)  Ball,  M33  .50  Caliber 
Projectiles 


All  Dinensions  in  Calibers 
□ne  Caliber  =  12.95nn 

Figure  2.  Schematic  of  .50  Caliber  Ball,  M33  Projectile 
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Figure  3.  Schematic  of  ,50  Caliber  Limited  Range  Training  Ammunition  (LRTA)  Design 


Figure  4.  Computational  Surface  Geometry  Model  of  LRTA 
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Figure  5.  Computed  and  Measured  Pitching  Moment  Coefficient  Versus  Mach  Number, 
M8  and  M33  Projectiles 
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Figure  6.  Computed  and  Meaeured  Normal  Force  Center  of  Pressure  Versus  Mach  Number, 
M8  and  M33  Projectiles 
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Figure  7.  Computed  and  Measured  Normal  Force  Coefficient  Versus  Mach  Number,  M8 
and  M33  Projectiles 
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Figure  8.  Computed  and  Measured  Magnus  Moment  Coefficient  Versus  Mach  Number,  M8 
and  M33  Projectiles 
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Figure  9,  Computed  Magnus  Force  Coefficient  Versus  Mach  Number,  M8  and  M33  Projec 
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Figure  10.  Computed  Magnus  Force  Center  of  Pressure  Versus  Mach  Number,  M8  and 
M33  Projectiles 
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Figure  11.  Computed  and  Measured  Pitch  Damping  Moment  Coefficient  Versus  Mach 
Number,  M8  and  M33  Projectiles 
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Figure  12.  Computed  Roll  Damping  Coefficient  Versus  Mach  Number,  M8  and  M33  Pro¬ 
jectiles 


20 


Figure  13.  Computed  and  Measured  Zero- Yaw  Drag  Coefficient  Versus  Mach  Number,  M8 
and  M33  Projectiles 


Figure  14.  Computed  Normal  Force  Coefficient  Versus  Body  Length,  LRTA  Configurations, 
Mach  2.6 
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Figure  15.  Computed  Normal  Force  Coefficient  Versus  Body  Length,  LRTA  Configurations, 
Mach  1.96 


Figure  16.  Computed  Normal  Force  Coefficient  Versus  Body  Length,  LRTA  Configurations, 
Mach  1.5 
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Figure  17.  Computed  Normal  Force  Center  of  Pressure  Versus  Body  Length,  LRTA  Con¬ 
figurations,  Mach  2.6 
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Figure  18.  Computed  Normal  Force  Center  of  Pressure  Versus  Body  Length,  LRTA  Con¬ 
figurations,  Mach  1.96 


23 


x/d 


Figure  19.  Computed  Normal  Force  Center  of  Pressure  Versus  Body  Length,  LRTA  Con¬ 
figurations,  Mach  1.5 
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Figure  20.  Computed  Roll  Damping  Coefficient  Versus  Body  Length,  LRTA  Configurations, 
Mach  2.6 
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Figure  21.  Computed  Roll  Damping  Coefficient  Versus  Body  Length,  LRTA  Configurations, 
Mach  1.96 
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Figure  22.  Computed  Roll  Damping  Coefficient  Versus  Body  Length,  LRTA  Configurations, 
Mach  1.5 
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Figure  23.  Computed  Roll  Damping  Coefficient  for  2-Fin  and  4-Fin  Configurations  Versus 
Body  Length,  Mach=1.96 


Figure  24.  Computed  Pitch  Damping  Moment  Coefficient  Versus  Mach  Number,  LRTA 
Configurations 
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Figure  25.  Computed  Zero- Yaw  Drag  Coefficient  Versus  Mach  Number,  LRTA  Configura¬ 
tion  #1 
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Figure  26.  Computed  Zero- Yaw  Drag  Coefficient  Versus  Mach  Number,  LRTA  Configura¬ 
tion  #2 
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Figure  27.  Computed  In-Flight  Mach  Number  Versus  Range,  LRTA  Configurations 
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Figure  28.  Computed  In-Flight  Gyroscopic  Stability  Factor  Versus  Range,  LRTA  Config¬ 
urations 
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Figure  29.  Computed  In-Flight  Nondimensional  Spin  Rate  Versus  Range,  LRTA  Configu 


rations 
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Figure  30.  Computed  In-Flight  Total  Angle  of  Attack  Versus  Range,  LRTA  Configurations 
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Figure  31.  Doppler  Radar  Measurement  of  In-Flight  Total  Drag  Versus  Range,  M33  and 
LRTA  Projectiles 


Figure  32.  Computed  Effect  of  Magnus  Moment  Coefficient  on  Dynamic  Stability,  LRTA 
Configuration 
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LIST  OF  SYMBOLS 


A 


reference  axea  of  projectile,  evaluated  as  7r(f^/4 


Cd 

Cdo 

CG 

Cl 

Cl. 

Cl 

Cl. 

Cm. 

Cm.. 

Cm,  +  Cm,, 
Cm 
Cn 
Cn. 

Cn.. 

Cn 

CPn 

CPm 

Cs 

d 

F 

A  A  A 

Es,F,G 

He 


^yy 

I 

M 

p 

p 

9 

Q 

T 

A 

s 


total  drag  coefficient 
zero-yaw  drag  coefficient 

center  of  gravity,  in  calibers,  measured  along  x-axis 
lift  coefficient 

lift  coefficient,  evaluated  as  ClIol,  a  in  radians 
roll  moment  coefficient 

roll  damping  coefficient,  evaluated  as  Cif{pdlVoo) 

pitching  moment  coefficient,  evaluated  as  Cm  lot 

Magnus  moment  coefficient,  evaluated  as  Cn/(pd/14o)/Q; 

pitch  damping  moment  coefficient 

pitching  moment  coefficient 

normal  force  coefficient 

normal  force  coefficient,  evaluated  as  Cn/oi 

Magnus  force  coefficient,  evaluated  as  C's/(pd/Yoo)a 

Magnus  moment  coefficient 

normal  force  center  of  pressure,  in  calibers,  measured  along  x-axis 

Magnus  force  center  of  pressure 

Magnus  force  coefficient 

reference  diameter  of  projectile 

aerodynamic  force 

inviscid  flux  vectors  of  transformed  gas  dynamic  equations 

inviscid  source  term  for  cylindrical  coordinate  formulation 

of  transformed  gas  dynamic  equations 

radial  moment  of  inertia 

transverse  moment  of  inertia 

reference  length 

aerodynamic  moment 

pressure 

projectile  spin  rate  (rad/s) 
dynamic  pressure,  pv^ 

vector  of  dependent  variables  of  gas  dynamic  equations 
distance  from  x-axis 

viscosity  vector  of  transformed  gas  dynamic  equations 
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Sc  viscous  source  term  for  cylindrical  coordinate  formulation 

of  transformed  gas  dynamic  equations 
Sd  dynamic  stability  factor 

Sd*  stability  coefficient 

Sg  gyroscopic  stability  factor 

V  velocity 

u,v,w  velocity  components  in  x,y^z  directions 

x,y,z  physical  Cartesian  coordinates 


Greek 

a 

£ 


P 


Symbols 

total  angle  of  attack;  pitch  angle;  yaw  angle 
total  energy  per  unit  volume  of  fluid 
density 

transformed  coordinates 


Subscript 

oo  free  stream  condition 

B  projectile  base  condition 


Note: 

Aerodynamic  force  coefficients  are  defined  as  \  ^ 
Aerodynamic  moment  coefficients  are  defined  as 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  AKL-TR-662 _ DateofReport - December  1994 - 

2.  Date  Report  Received. - — - - - 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ _ _ — - 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) - - - - - 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. - - - 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.)  - - - - 


Organization 


CURRENT  Name 

ADDRESS  _ _ _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ - 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


Department  of  the  Army 
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BUSINESS  REPLY  MAIL 

FIRST  CUSS  PERMIT  NO  0001,  APG,  MD 

Postage  will  be  paid  by  addressee 

Director 

U.S.  Army  Research  Laboratory 

ATTN:  AMSRL-OP-AP-L 

Aberdeen  Proving  Ground,  MD  21005-5066 
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